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Abstract. A formalism for study of spectral correlations in non-Gaussian, 
unitary invariant ensembles of large random matrices with strong level con- 
finement is reviewed. It is based on the Shohat method in the theory of 
orthogonal polynomials. The approach presented is equally suitable for de- 
scription of both local and global spectral characteristics, thereby providing 
an overall look at the phenomenon of spectral universality in Random Ma- 
trix Theory. 



1. Introduction: Motivation and Basic Results 

1.1. UBIQUITY OF INVARIANT RANDOM MATRIX MODELS 

Random matrices [1, 2] have been introduced in a physical context since 
the pioneering works by Wigner [3] and Dyson [4, 5]. Initially proposed as 
an effective phenomenological model for description of the higher excita- 
tions in nuclei [6] they found numerous applications in very diverse fields 
of physics such as two dimensional quantum gravity [7], quantum chro- 
modynamics [8], quantum chaos [9], and mesoscopic physics [10, 11]. One 
can state that from the standpoint of the mathematical formalism all these 
fields are pooled by the Random Matrix Theory (RMT). Such an ubiquity 
of random matrices owes its origin to the exclusive role played by symme- 
try. The most amazing evidence to this fact comes from invariant random 
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matrix models which will be the focus of this review. The main feature of 
these models lies in that they discard (irrelevant) microscopic details of 
the physical system in question, but they do properly take into account 
its underlying fundamental symmetries. In accordance with the very idea 
of the construction of invariant matrix models, they do not relate to any 
dynamical properties of the physical object under study: general symmetry 
requirements alone lead to appearance of knowledge about the system. As 
far as the symmetry constraints follow from the first principles, the Random 
Matrix Theory turns out to be a general and powerful field-theoretical ap- 
proach leading to a unified mathematical description of the quite different 
physical problems mentioned above. 

A great variety of invariant random matrix ensembles can be assigned 
to three irreducible symmetry classes [5]. To specify them, we consider the 
typical line of arguing used in applications of the Random Matrix Theory 
to the disordered quantum mechanical systems, where it was first invented 
by Gor'kov and Eliashberg [12]. Since in this situation the microscopic 
Hamiltonian 7i is rather intricate, the integration of exact equations is 
impossible. It is therefore useful to appeal to statistical description by con- 
jecturing that the operator Ti can be modelled by an N x A?" random matrix 
H whose eigenvalues and eigenvectors reproduce statistically the eigenlevels 
and eigenfunctions of the real microscopic Hamiltonian in the thermody- 
namic limit which corresponds to the matrix dimension N going to infinity. 
With this conjecture accepted, an ensemble of large random matrices H, 
characterized by the joint distribution function P [H] of the matrix ele- 
ments Hy of the corresponding Hamiltonian 7i, becomes the main object 
of study. Once the primary role of symmetry is postulated in the RMT- 
approach, the matrix H must adequately reflect the symmetry properties of 
the physical system under study. The matrix H is chosen to be real symmet- 
ric if the underlying physical system possesses time-reversal and rotational 
invariance. Systems with broken time-reversal symmetry are characterized 
by a Hermitean matrix H, while systems with conserved time-reversal sym- 
metry but with broken rotational invariance are described by a self-dual 
Hermitean matrix. These three symmetry classes referred to as the orthog- 
onal, unitary and symplectic symmetry classes, respectively, can be char- 
acterized by a symmetry parameter (3 equals the number of independent 
elements in the off-diagonal entries of the matrix H. The parameter (3 = 1 
for a real symmetric matrix (orthogonal symmetry), (3 = 2 for a Hermitean 
matrix (unitary symmetry), and (3 = 4 for a self-dual Hermitean matrix 
(symplectic symmetry). 

So far we did not yet specify the form of the joint distribution func- 
tion P [H] of the matrix elements Hjj. By definition, the invariant random 
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matrix ensembles are characterized by 

P[H] = ^Le X p{-/3TrV[H]}. (1) 

Here the function V [H] should ensure the existence of the partition func- 
tion Zn, which is defined by the normalization condition / P [H] d [H] = 1 
with the elementary volume d [H] depending on the symmetry of the matrix 
H. For (3 = 1 the volume element d [H] = Yli<jdHij, for (3 = 2 the vol- 
ume element d [H] = Yii<j dReHij IL<j dlmtlij, while for (3 = 4 it equals 

d [H] = Ui<j dn\f Ul=i Ui<j dUff- The presence of the trace in Eq. (1) 
leads to the invariance of the probability density P [H] d [H] under the simi- 
larity transformation H — ► TZ^HTZp with TZp being an orthogonal, unitary 
or symplectic N x N matrix for (3 = 1, 2 or 4, respectively. In turns, the 
invariance built in the probability density P [H] d [H] implies that there is 
no preferential basis in the space of matrix elements. From the physical 
point of view, this means that invariant matrix models given by Eq. (1) 
are applicable to particular regimes of a physical system where (i) all the 
normalized linear combinations of the eigenstates have similar properties, 
and where (ii) the dimensionality is irrelevant. In disordered systems this 
is just a metallic state where the typical electron states are extended and 
hence structureless. 

Notice that up to this point we have no constraints allowing us to 
uniquely choose the function V [H] (referred to as "confinement poten- 
tial"). However, if we impose the additional requirement that the entries 
of the random matrix H be statistically independent of each other, we 
immediately arrive at the Gaussian Orthogonal, Unitary and Symplectic 
Ensembles (GOE, GUE, GSE) which are characterized by the quadratic 
confinement potential V [H] oc H 2 . This particular form of confinement po- 
tential leads to significant mathematical simplifications which allowed the 
complete treatment of these three ensembles many years ago [13]. 

It is remarkable that (even for non-Gaussian distributions of P [H] ) the 
invariant random matrix ensembles possess a great degree of mathemat- 
ical tractability, and, what is more important, they have a high physical 
relevance, being much more than just a mathematical construction. In this 
respect we mention that in the physics of disordered systems the applicabil- 
ity of Gaussian invariant random matrix ensembles to description of weakly 
disordered systems has been proven by Efetov [14] by using the nonlinear 
c-model. In that study [14] the statistical properties of energy levels for 
metallic particles with volume imperfections were considered by solving the 
Schrodinger equation with nonperturbative averaging over the random po- 
tential configurations within the framework of the supersymmetry method. 
Random Matrix Theory appears there as a zero-dimensional version of a 
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more general microscopic nonlinear a-model, thereby proving the validity 
of the basic principles used for an RMT phenomenological description of 
energy levels of noninteracting electrons confined in a restricted volume. 
This connection takes place at times much larger than the ergodic time 
needed for diffusive particle to completely and homogeneously fill the avail- 
able volume of the sample provided it is in the metallic regime characterized 
by a dimensionless Thouless conductance g> 1. 

Let us point out that the choice of quadratic confinement potential 
V [H] can hardly be justified. Indeed, it was understood from the very 
beginning [13] that the requirement of statistical independence of the ma- 
trix elements Hjj is not motivated by the first principles and, therefore, 
the important problem of elucidating the influence of a particular form of 
confinement potential on the predictions of the Random Matrix Theory 
developed for Gaussian Invariant Ensembles had been posed already in the 
sixties. Despite this fact, considerable progress in study of spectral proper- 
ties of non-Gaussian random matrix ensembles was achieved almost thirty 
years later when RMT experienced a great renaissance due to new ideas in 
the physics of disordered/chaotic systems which had led the birth of meso- 
scopic physics, as well as due to a penetration of Random Matrix Theory 
to quantum chromo dynamics (QCD). In the latter field, the Random Ma- 
trix Theory turned out to be a useful tool for understanding the spectral 
properties of low-lying eigenvalues of the Dirac operator. The idea of in- 
troducing the RMT-approach in QCD is very similar to that in the physics 
of disordered systems and is based on the conjecture [8, 15] that the spec- 
tral density of the Dirac operator very close to the spectrum origin should 
depend only on the symmetries in question. One startling consequence of 
this conjecture is that the spectral density of the Dirac operator near the 
origin need not be computed within the framework of the gauge theories at 
all, but it can be extracted from much simpler random matrix ensembles 
reflecting the symmetry of the problem. The matrix models appearing in 
the context of QCD are manifestly non-Gaussian possessing an additional 
chiral structure 1 [16]. Another motivation for studying the spectral prop- 
erties of non-Gaussian random matrix ensembles comes from the theories 
of 2D quantum gravity [7]. 

1.2. NON-GAUSSIAN RANDOM MATRIX ENSEMBLES AND 



throughout the paper we consider non-chiral non-Gaussian random matrix ensem- 
bles. It can be shown that an arbitrary chiral matrix model can be reduced to an auxiliary 
one without chirality; see, for instance, Ref. [17]. 
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PHENOMENON OF SPECTRAL UNIVERSALITY 

The examples above clearly demonstrate an important role the non-Gaussian 
random matrix ensembles play in different physics theories, and serve as a 
compelling evidence of the necessity to have a powerful method for study 
of their spectral properties which could be equally applicable to rather 
different probability measures P [H] . During (mostly) the recent decade a 
number of methods were developed in order to treat non-Gaussian random 
matrix ensembles. All of them can schematically be related to two groups. 

1.2.1. Global Universality 

The first group includes different approximate methods useful to explore 
global spectral characteristics which manifest themselves on the scale of 
n S> 1 eigenlevels. Among these methods there are (i) the mean-field 
approximation proposed by Dyson [4, 18] which allows one to compute 
the density of levels in random matrix ensemble; (ii) the Schwinger-Dyson 
loop equations' technique [19, 20, 21] that had led to discovery of the phe- 
nomenon of the global spectral universality; (iii) the method of functional 
derivative of Beenakker [22, 23], and (iv) the diagrammatic approach of 
Brezin and Zee [24] whose development enabled their authors to study the 
phenomenon of global universality [19] in more detail as well as to generalize 
it in the context of mesoscopic physics. 

It was found that contrary to the one-point spectral characteristics 
(such as one-point Green's function or level density) which essentially de- 
pend on the measure P [H], i.e. on the explicit form of the confinement po- 
tential, the functional form of (connected) two-point correlators becomes 
insensitive to the details of confinement potential upon smoothing over 
the scale which is much larger than the mean level spacing An but much 
smaller than the scale of the entire spectrum support. This is the essence of 
the phenomenon of global spectral universality. For example, the smoothed 
connected density-density correlator computed for the matrix model Eq. 
(1) is given by the universal function 



where A 7^ A'. This form of p c is valid for random matrices whose spec- 
trum is supported on a single, symmetric interval (—Vn,+T>n). Here, the 
universality implies that all the information about the particular form of 
V [H] is encoded into pi N ^ only through the end point of the spectrum 
support. 

The methods mentioned above, being applicable to study of the spectral 
correlations in the long-range regime for all three symmetry classes, leave 






6 



aside the fine structure of eigenvalue correlations manifested on the scale 
of the mean eigenvalue spacing. In this sense, these approaches are less 
informative as compared to the method of orthogonal polynomials [1] which, 
along with the supersymmetry approach [25, 26], enters the second group. 

1.2.2. Local Universality 

At present, the orthogonal polynomial technique, originally developed by 
Gaudin and Mehta [27], seems to be the most powerful one furnishing us 
the possibility to probe both the global spectral characteristics (which are 
of importance in computing the integral spectral properties) and the local 
spectral correlations (which describe a dynamics of underlying physical 
system) for arbitrary symmetry classes. 

(i) The early attempts [28, 29, 30] to go beyond the Gaussian distribu- 
tion of P [H] were concentrated on unitary invariant, U (N), matrix ensem- 
bles ((3 = 2) associated with classical orthogonal polynomials. It was found 
in Refs. [28, 29, 30, 31] that (in the thermodynamic limit N — ► oo) the 
scalar two-point kernel (in terms of which all n-point correlation functions 
are expressible, see Sec. 2.2) computed in the bulk of the spectrum, i.e. far 
from the end points of the eigenvalue support, follows the sine law 



-Kbulk (s, s') 



sin [-7T (s — s')] 

IT (s — s') 



(3) 



that inevitably leads to the famous Wigner-Dyson level statistics [32] inher- 
ent in nuclear physics, weakly disordered and chaotic systems. This form of 
the scalar kernel corresponds to the bulk scaling limit, when the initial spec- 
tral variables A, A' are measured in the units of the mean eigenvalue spacing 
Ajv, so that the scaled variable s = A/A^r- The same form of the scalar 
kernel Eq. (3) obtained for different ensembles associated with classical or- 
thogonal polynomials was the first evidence to the phenomenon currently 
known as the local spectral universality. Later, it was shown in Ref. [31] that 
spectral properties of orthogonal O (N) and symplectic Sp (./V) {(5 = 1 and 
(3 = 4) invariant ensembles associated with the weights of classical orthog- 
onal polynomials 2 follow the predictions of GOE and GSE, respectively. 

A further significant progress in the field came with the works [35, 36] 
whose authors considered spectral properties of U (N) invariant random 
matrix ensembles associated with strong symmetric confinement potentials 
of the form V [H] = H 2 + 7 H 4 (7 > 0) and V [H] = £*=i a fc H 2fc (a p > 0), 

2 In this situation the spectral correlations are expressible through the 2x2 matrix 
kernel, that can be computed by using so-called skew orthogonal polynomials [33]. In 
fact, as was recently shown in Ref. [34], the ensembles with (3 = 1 and f3 = 4 can be also 
treated without appealing to the skew polynomials. 
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respectively. Both works, based on different conjectures about the func- 
tional form of asymptotics of polynomials orthogonal with respect to a 
non-Gaussian measure, restricted their attention to the spectrum bulk, 
where the two-point kernel was shown, once again, to follow the sine law, 
Eq. (3). A rigorous treatment of a richer class of U (N) invariant random 
matrix ensembles related to the Freud and Erdos-type orthogonal polyno- 
mials was given in Refs. [37, 38]. It was proven there that the universal 
sine law for the two-point kernel holds for a wide class of monotonic (not 
necessarily of polynomial form) non-singular confinement potentials V (A) 
which increase at least as |A| at infinity, and can grow as or even faster than 
any polynomial at infinity. Confinement potentials satisfying these proper- 
ties are referred to as strong confinement potentials. This definition takes 
its origin in the limits [37] of spectral universality and is non-accidently 
connected to the problem of determinate and indeterminate moments [39] . 
Also, it was demonstrated in Ref. [38] that an intimate connection exists 
between the structure of Szego functional [40] entering the strong point- 
wise asymptotics of orthogonal polynomials and the mean-field equation 
by Dyson for mean level density that has been derived in Ref. [38] within 
the framework of orthogonal polynomial technique. 

(ii) All the random matrix ensembles treated in Refs. [28, 29, 30, 31, 
35, 36, 37, 38] were characterized by strong confinement potential with no 
singularities. However, (logarithmic) singularities do appear when one con- 
siders chiral matrix ensembles, arising in the context of QCD [8], in the the- 
ory of mesoscopic electron transport [41] and in description of electron level 
statistics in normalconducting-superconducting hybrid structures [42]. An 
example of a rather general (though non-chiral) random matrix ensemble 
possessing a log-singular level confinement is given by the distribution 

P[H] = ^HdetH| Q/3 exp{-/31W [H]}, (4) 

where the function V [H] is a well behaved function which has not singu- 
lar points. Ensemble Eq. (4), being a natural generalization of the matrix 
ensemble proposed by Bronk [30], was first considered in the situations 
where associated orthogonal polynomials were classical [30, 43, 44, 45]. For 
quadratic confinement potential and (3 = 2 one obtains that in the vicinity 
of the singularity, A = 0, the scalar kernel satisfies the Bessel law 

„ / A _ 7T , a 1/2 Ja+1/2 (7[g) J q -1/2 C^fO ~ Jg-1/2 C^f) Jg+1/2 (^0 
-ftorig [S, S ) — ^ [ ss ) s — s' ' 

" " ^ , (5) 

where s and s' are scaled by the level spacing Atv (0) near the spectrum 
origin, s = A/A^ (0), and a > —1/2. This scaling procedure is referred to 
as the origin scaling limit. Extensions to two other symmetry classes, as well 
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as to the chiral matrix ensembles, can be found in Refs. [43, 44, 45, 46]. 
An important breakthrough in understanding the universal character of 
this kernel was given in Refs. [47, 17]. These authors, guided by QCD 
applications, have shown that the Bessel kernel Eq. (5) is again universal, 
being independent of the details of strong confinement potential V [H]. An 
alternative proof of universality that holds more generally was presented in 
Ref. [48]. 

(iii) The third type of universal correlations takes place near the soft 
edge of the spectrum support, which is of special interest in the models 
of two-dimensional quantum gravity. The first study of the level density 
near the end point of the spectrum support is due to Wigner [49]. More 
comprehensive description of the tails of the density of states was done in 
Ref. [50] . It was shown there that at (5 = 2 a universal crossover occurs from 
a nonzero density of states to a vanishing one that is independent of the 
confining potential in the soft-edge scaling limit. Later it was demonstrated 
[51] that eigenvalue correlations in the U (N) invariant matrix ensembles 
with quartic and sextic confinement potentials are determined by the scalar 
kernel obeying the Airy law [52] 

Ai (a) Ai' (gQ - Ai (s') Ai' (s) 

-Ksoft (s, s ) = , (bj 

s — s' 

suggesting that the Airy kernel should be universal as well. (Here the rescal- 
ing s oc iV 2 / 3 (A/Pjy — 1) determines the soft -edge scaling limit). This con- 
jecture has been proven in Ref. [53], where it was also shown that the Airy 
correlations, being universal for a class of matrix models with monotonic 
confinement potential or with that having light local extrema, are indeed 
a particular case of more general universal multicritical correlations [53]. 

1.2.3. Universality in a Broader Context 

Are the universal scalar kernels in unitary invariant random matrix models 
with strong level confinement exhausted by the universal sine, Bessel and 
Airy laws given by Eqs. (3), (5) and (6) above? The present state of the art 
leads us to the negative answer. Indeed, by adding a singular component to 
the strong level confinement we may either accumulate a finite number of 
eigenvalues near the singular point A s i ng or repel them from it. Such a rear- 
rangement of eigenlevels will lead to emergence of a new scalar two-point 
kernel in the vicinity of the singular point of the spectrum. (In particular 
case of the logarithmic singularity this kernel will follow the Bessel law 
Eq. (5)). Generically, it is natural to expect that the functional form of 
the kernel near A s i ng will be sensitive to the particular type of the singu- 
lar deformation. However, the new two-point kernel will be insensitive to 
the details of the background component of the confinement potential, as 
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it takes place in the case of the Bessel kernel. In this sense, one can still 
say about a phenomenon of universality. One of the latest evidences to this 
fact can be found in Ref. [54] where deformed ensembles of large random 
matrices associated with massive Dirac operators were considered. 

By the same token, the universal global spectral correlator expressed 
by Eq. (2) is not the only one arising in Random Matrix Theory. As it was 
already stressed, the universal function Eq. (2) is inherent in matrix mod- 
els with spectra possessing a single connected support. Correspondingly, 
ensembles of large random matrices with eigenvalue densities having more 
than one-cut support will give rise to the novel global spectral correlators 
whose universality classes (for a given symmetry parameter (3) are charac- 
terized entirely by the number of cuts in the support of spectral density. 
This was explicitly demonstrated in the recent studies [20, 55] by means of 
the loop equation technique. 

1.3. THE AIM 

At this point, it is appropriate to notice that all the mentioned above (uni- 
versal) results for both global and local eigenvalue correlations have been 
obtained by using different methods, each of them was only suitable for a 
particular problem under consideration. Any deformations of P [H] (which 
will preserve its invariance) would cause principal difficulties in elucidating 
the influence of these deformations on spectral properties of corresponding 
random matrix ensembles. The goal of this paper is to represent a general 
method (recently developed in a series of publications [53, 48, 56]), which is 
equally suitable for study of both local and global eigenvalue correlations, 
and easily leads to generalizations. The approach we introduce is based on a 
simple and elegant idea [57] by J. Shohat (which goes back to 1930), provid- 
ing a detailed description of the spectral properties of non-Gaussian U (N) 
invariant random matrix ensembles through the analysis of the three-term 
recurrence equation for associated orthogonal polynomials. We show that 
for the most situations of interest, the knowledge of the large-A behav- 
ior of the coefficients in the recurrence equation is sufficient to directly 
reconstruct the local eigenvalue correlations of arbitrary order, as well as 
to explore the global spectral statistics. In the case of a non-singular, well 
behaved confinement potential, the knowledge of such a large-A behav- 
ior of the recurrence coefficients is equivalent, in fact, to a knowledge of 
the Dyson density of states for the corresponding random matrix ensem- 
ble. The latter assertion leads to a rather unexpected conclusion: Once the 
Dyson density of states (which is a rather crude one-point spectral char- 
acteristics) is available, the scalar kernel (and hence the n-point spectral 
correlators) can immediately be recovered through the solution of a certain 
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second-order differential equation (See Sec. 4.3). 

We believe that this method offers not only new computational po- 
tentialities, but also provides a different, overall look at the problem of 
eigenvalue correlations in unitary invariant random matrix ensembles in 
arbitrary spectrum range and in arbitrary scaling limits. It seems that to- 
gether with the very recent works [58, 59] establishing a precise connection 
of the scalar kernel for random matrix ensembles with U (N) symmetry 
with the 2x2 matrix kernels in ensembles with O (N) and Sp (JV) symme- 
tries, the formalism to be reviewed below gives a rather complete solution 
of the problem of eigenvalue correlations in invariant matrix models with 
strong level confinement. 

The review is organized as follows. Section 2 contains a brief description 
of the Gaudin-Mehta calculational scheme, that introduces the orthogonal 
polynomials as a tool for exact evaluation of n-point correlation functions in 
U (iV) invariant random matrix ensembles. In Section 3 the Shohat method 
in the theory of orthogonal polynomials is presented. Section 4 is devoted 
to a detailed study of spectral properties of large Hermitean random matri- 
ces with a single connected eigenvalue support. In Section 5 we extend this 
analysis to random matrices with eigenvalue gap. Section 6 contains con- 
clusions. The most lengthy calculations are collected in three Appendices. 

2. Elements of the Gaudin Mehta formalism 

2.1. INVARIANT RANDOM MATRIX MODEL IN EIGENVALUE 
REPRESENTATION: TWO INTERPRETATIONS 

The invariance of the distribution function P [H] implies that different ma- 
trices with the same eigenvalues have the same probability of occurring. 
To study spectral characteristics of an invariant random matrix model it is 
convenient to integrate out "auxiliary" angular variables in the construction 
P [H] d [H] in order to get the matrix model in the eigenvalue representa- 
tion. To proceed with this, we have to pass from the integration over inde- 
pendent elements Hy of the matrix H to the integration over the smaller 
space of its N eigenvalues {A}, calculating the corresponding Jacobian J. 

Let us introduce the matrix Tip that diagonalizes the random matrix 
H, so that H = TV^ATlp and A = diag (Ai, Ajv), and consider the 
infinitesimal variation of H, 

<5H =TT^ (llpSTl^A + 5A + ASTlpTl^ lip = K^ 1 (5A - i [A, 5s])Tlp. 

2 (7) 

Here we have denoted 5s = iSHpUp 1 . The norm, ||«5H|r = Tr(<5H) 2 , of the 
infinitesimal variation of H is 

||<5H|| 2 = Tr (5 A) 2 - 2i Tr ( [5 A, A] 5s) + 2 Tr (-<5sA(5sA + (5s) 2 A 2 ) . (8) 
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The second term in the last expression vanishes as the matrices <5A and A 
are diagonal. We then find 

||5H|| 2 = (**0 2 + £ ( A * - x j) 2 \ Ss a\ 2 ■ ( 9 ) 

The independent variables are the variations of the eigenvalues <5Aj and Ssij 

for 0=1, ReSsij, ImSsij for = 2, or 8s\f with a = 0,1,2,3 for = 4, 
i < j. Then, from Eq. (9) we obtain the Jacobian J of the transformation 
H ^ {\i,TZfj} that is equal to V det G, where G is the metric tensor with 
det G = Y\i-Lj |Aj — \j\ . Hence 

J = \{\\-\ 3 f = \^{\)f (10) 

i<j 

with A (A) being the Vandermonde determinant. The construction Eq. (10) 
is also known as the Jastrow factor. 

Combining Eqs. (1) and (10) we arrive at the famous expression for the 
joint probability density function of the eigenvalues {A} of the matrix H: 

N 

p(X 1 ,...,X N ) = Z^exp{-0Y / V(X l )}\A(X)f (11) 

i=i 

= Z^e W {-P[J2 V ( X i)-J2 ln \^- X j\]}- ( 12 ) 

i i<j 

The probability distribution given by Eq. (12) has the form of a Gibbs dis- 
tribution for a classical one—dimensional system of N "particles" at "posi- 
tions" Aj confined by the external one-body potential V (A) and interacting 
with each other through the pairwise logarithmic law originated from the 
Jacobian J of the transformation H ^{Xi,lZp}. The symmetry parame- 
ter plays the role of the equilibrium temperature. Such an interpretation 
of Eq. (12) gives rise to approximate mean-field methods in the Random 
Matrix Theory. 

For invariant matrix ensembles with unitary symmetry {0 = 2), the 
probability distribution in the form of Eq. (11) can alternatively be related 
to a system of fictitious non-interacting fermions, that can be described 
by effective one-particle Schrodinger equation [53]. This equation is a cor- 
nerstone of the method under review. Although such a simple, transparent 
interpretation cannot be ascribed to p (Ai, Xn) for two other symmetry 
classes, = 1 and = 4, a recently discovered deep connection [58, 59] 
between orthogonal, unitary and symplectic ensembles makes the unitary 
invariant ensembles of random matrices to be the central and most impor- 
tant. 
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2.2. ORTHOGONAL POLYNOMIALS' TECHNIQUE: (3 = 2 

For (3 = 2 the structure of Eq. (11) enables us to exactly represent all 
the statistical characteristics of the spectrum in the terms of polynomials 
orthogonal with respect to a non-Gaussian measure da. To demonstrate 
this, it is useful to write down the joint probability density function of the 
eigenvalues {A} in the form 



p(Ai,...,Ajv) = *g(Ai,...,A JV ), (13) 

N 

* (Ai,...,A JV ) = Z" 1/2 exp{-^y(A i )}A(A). (14) 



=i 



One can see that can be thought of as a wave function of N fictitious 
non-interacting fermions. Namely, noting that 

A(A)= J] (A i -A j ) = det||^'- 1 |, (15) 

i>j=l 

and taking the linear combinations of the columns of the initial matrix with 
entries X{~ 1 , one can reduce this matrix to the matrix whose entries are 
arbitrary polynomials Pj-i (Aj) of degrees j — 1 = 0,1,..., A/" — 1, 

A (A) = det \\Pj-i (Aj)|| . (16) 

Further, choosing these polynomials to be orthogonal with respect to the 
measure da (A) = exp {— 2V (A)} d\, 



da(X)P n (X)P m (X) = S nm , (17) 

Ae(— oo,+oo) 

we arrive at the conclusion that the function m Eq. (14) can be repre- 
sented as a Slater determinant 

*o(Ai,...,Aiv) = det ||^_i (Ai)|| (18) 

that formally describes the system of ./V fictitious non-interacting fermions 
located at "spatial" points Aj and characterized by the set of orthonormal 
"eigenfunctions" 

<p n (X) =P n (X)exp{-V(X)}, (19) 

/+oo 
dXifn (A) ip m (A) = 5 nm . (20) 
-oo 
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Bearing in mind the representation Eq. (13), and taking into account Eqs. 
(17) - (19), we readily get that 

p (Ax, . . . , X N ) = ^ det \\K N (A,, . . =1 _ N , (21) 

where 

Af-l 

K N {X,X') = ^(A)^(A') (22) 

is the scalar two-point kernel. Making use of the Christoffel-Darboux the- 
orem [60] (see Eq. (32) below), the two-point kernel can be reduced to the 
form 

r, a wn yjy (AQ i£N-i (A) - ip N (A) fN-i (AQ 
K N {X,X)=c N _ , (23) 

convenient for the further analysis. In Eq. (23) cjy is the coefficient in 
the three-term recurrence equation (see Eq. (28) below) for the set P n of 
the polynomials orthogonal with respect to the measure da. This formula 
simplifies significantly RMT calculations, since the "eigenfunctions" (fN 
with large "quantum numbers" N » 1 entering Eq. (23) can be replaced 
by their large-iV asymptotics. 

The n-point correlation function is determined in the Random Matrix 
Theory by the formula [1] 

TV! it r +oc 
R n (X 1 ,...,X n )= [] / dXkp (Ai, . . . , Aat) . (24) 

It describes the probability to find n levels around each of the points 
Ai, . . . , A n when the positions of the remaining levels are unobserved. The 
multiple integrals in the last equation can exactly be calculated by using 
the representation Eq. (21). The result of the integration reads [1] 

R n {X 1 ,...,X n ) = det\\K N (X i ,X j )\\ ij=1 n . (25) 

Equation (25) implies that the knowledge of the scalar two-point kernel 
Kn (A, A') is sufficient to calculate the energy level correlation function of 
any order. In particular, the averaged density of states is expressed as 

(u N (X)) = R 1 (X) = K N (X,X). (26) 

Analogously, Eq. (25) leads to the following expression for connected "density- 
density" correlation function pi N ^ = (i/jy (A) (A')) — (vn (A)) (vn (A'))> 



pf) = S (A - A') R 1 (A) + R 2 (A, A') 

= 5(A-A')^(A,A)-^(A,A'). 



(27) 
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Thus, all the nontrivial information about eigenlevel correlations is con- 
tained in the squared two-point kernel — (A, A'). 

3. The Shohat Method: General Relations 

Equations (13) and (18) suggest that the joint distribution function of N 
eigenvalues of a U (N) invariant random matrix ensemble can be interpreted 
as a probability of finding N fictitious non-interacting fermions to be con- 
fined in a one-dimensional space. The effective one-particle Schrodinger 
equation for the wave functions ip n (A), Eq. (19), of these fictitious fermions 
can be derived by mapping a three-term recurrence equation for orthogo- 
nal polynomials, Eq. (17), onto a second-order differential equation. The 
method of reducing a recurrence equation to a differential equation is es- 
sentially due to J. Shohat who proved in 1939 that orthogonal polynomials 
associated with exponential weights satisfy a second-order differential equa- 
tion [57]. Much later the Shohat method has got a further development in 
the work by Bonan and Clark [61]. By now, rather extended mathematical 
literature exists on this subject [62, 63, 64]. 

Let as consider a set of polynomials P n (A) orthogonal on the entire real 
axis with respect to the measure da (A) = exp {— 2V (A)} dX. If V (A) is an 
even function 3 , V (—A) = V (A), this set of orthogonal polynomials can be 
defined by the recurrence equation 

AP n _i (A) = c n P n (A) + c„_iP n _ 2 (A) , (28) 

where the coefficients uniquely determined by the measure da. 

In order to derive the differential equation for the wave functions (p n (A) = 
P n (A) exp {— V (A)}, we note that the following identity takes place, 

= A n (A) P n _! (A) - B n (A) P n (A) , (29) 

with functions A n (A) and B n (A) to be determined from the following con- 
sideration. Since dP n (A) /dX is a polynomial of the degree n — 1, it can be 
represented [60] through the Fourier expansion in the terms of the kernel 

Qn (t, X) = Efc=d p k (A) Pk (t) as 

Integrating by parts in the last equation we get that 

dPn(X) f , , x ^ , fdV dV\ „ . , 
-±> = 2j l la( t )Q nit .X ) (---)p n (t). (31) 

3 For asymmetric confinement potentials the recurrence equation takes the form 
(A — b n ) Pn-i (A) = c„P„ (A) + c n _iP„_2 (A) . The additional parameter b n can easily 
be incorporated into the calculational scheme. 
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Now, making use of the Christoffel-Darboux identity [60] 
n-l 

P, (\) P, (t) = r- {Z> { ' 

t-X 



Qn (t, A) = £ ft (A) ft (t) = c ^O^-lW-fnW^-lO , 



(32) 



fc=0 

we are led to 



dA " y w t - A 

x [P n (t) P n _! (A) - P n (A) P n _! (*)] . (33) 

Comparison of this expression with Eq. (29) yields 

A n {\) = 2c n J da (t) V ' {t ]Z V x ' {X) P*(t), (34) 

B n (X) = 2c n J da (t) V ' (t | ~ ^ (A) P n (t) P n _x (t) . (35) 

Now one can obtain the exact differential equation for the eigenfunctions 
<p n . Differentiating Eq. (29), consequently applying Eqs. (29) and (28), and 
taking into account Eq. (19) we derive after somewhat lengthy calculations 



<P<Pn (A) m difn{X) 

— ^2 J-n (A) — — V y n (A) <p n (A) = 0, (36) 



where 



*•<*> ^ it < 37 » 

„ / * \ dB„ c„ . . „ / „ dV 1 dA n 

Qn (A) = -rf + ^-A n A n -i-B n [B n + 2— + 
dX Cn-i \ dX 



d 2 V fdV\ 2 1 dA n dV 



A n dX 

+ ~dX 2 ~~{~dx) ~A^~dX~dX' ^ 
When deriving Eqs. (36), (37) and (38) we made use of the sum rule 

B n + B n _i - — A»-i = -2^, (39) 

Cn-l OA 

directly following from Eqs. (34), (35), (28) and from oddness of dV/dX. 

Equation (36) is valid for arbitrary n. Previously, an equation of this 
type was known in the context of the Random Matrix Theory only for 
GUE, where V (A) = A 2 /2. For such a confinement potential both functions 
A n (A) and B n (A) can readily be obtained from Eqs. (34) and (35), and are 
given by A n (A) = 2c n and B n (A) = 0. Taking into account that for GUE 
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c n = (n/2) 1/2 , we end up with T n (A) = and Q n (A) = 2n + 1 - A 2 . This 
allows us to interpret ip n (A) as a wave function of the fermion confined by 
a parabolic potential, 

\ 2 (A) + (2n + 1 - A 2 ) ^ UE (A) = 0. (40) 

The effective Schrodinger equation (36) applies to general non-Gaussian 
random matrix ensembles as well, although the explicit calculation of T n (A) 
and Q n (A) in this situation is a rather complicated task. In two cases of 
relatively simple measures with V (A) = A 4 /8 + g3A 3 /6 + g2A 2 /4+giA/2 and 
V (A) = A 6 /12 the functions T n (A) and Q n (A) are known explicitly [62]. 
Significant simplifications, however, arise in the limit n = N S> 1, which is 
just a thermodynamic limit of the Random Matrix Theory representing for 
us the most interest. 



4. Random Matrices with Single Eigenvalue Support 

In this Section we will be interested in the study of eigenvalue statistics for 
unitary invariant non-Gaussian large random matrices characterized by a 
distribution function P [H] given by Eq. (4). The confinement potential 
associated with this model is 

Va (A) = v (A) — a log | A| . (41) 

Here v (A) is the regular part of level confinement 

^) = £t^ (42) 

k=l 

with a positive leading coefficient, d p > 0; the signs of the rest of the can 
be arbitrary but they should lead to an eigenvalue density supported on a 
single connected interval, {A} € (— T>n,+Vn). The parameter a > —1/2 is 
the strength of the logarithmic singularity. 

In Subsection 4.1 below, we demonstrate how the one-point spectral 
characteristics (density of states) can be obtained by making use of the 
recurrence equation (28). In Subsection 4.2, we turn to the study of the 
smoothed connected "density-density" correlator, also starting with recur- 
rence equation (28). Finally, in Subsection 4.3, we obtain the universal 
scalar kernels in the origin, bulk and soft-edge scaling limits by solving the 
effective Schrodinger equation for fictitious fermions. 



4.1. MACROSCOPIC LEVEL DENSITY FROM RECURRENCE EQUATION 

We start with an explanation of the main idea of the derivation to make 
clear all the subsequent calculations. Here we mainly follow Ref. [65]. Our 
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basic observation is that in the large-iV limit the density of states 4 vn (A) 
consists of two parts, 

MA) = z4 m °° th (A) + ^ c (A). (43) 

The smooth part i/ s ™ oth - (A) contributes to different integral characteristics 
determined by the density of states, while the oscillating part does not, 
because any integration will level the oscillating features. Then, for some 
smooth, well behaved, even function / (A) we have 

/ + °° dXf (A) VN (A) N =?° I dXf (A) ^ m °° th (A) • (44) 

J — oo JAGsupport 

Let us implement this scheme by choosing (without any loss of general- 
ity) / (A) = A 2s , with s being a positive integer. This choice is possible due 
to the evenness of vn (A). By definition Eq. (26), we obtain from Eq. (23) 

v N (A) = c N exp {-2V a (A)} (p£>' (A) P^ (A) - pg (A) p£> (A)) . 

(45) 

Having in mind the relation Eq. (29), the sum rule Eq. (39) and the defi- 
nition Eq. (19) we come down to 



MA) = °n (A) (y^-i (A)) + ^4 Q li (A) (# (A)) (46) 

-V$ (A) <P%U (A) (— (A) + BP (A) - BW (A))" . 

Here the upper index (a) is used to reflect the presence of a log-singular 
component in confinement potential V a (A), Eq. (41). Thus, the level density 

is expressed through the wave functions tpff (A), and through the functions 

A { p (A) and Bp ] (A), given by Eqs. (34) and (35). 

For further convenience we introduce two quantities, 

A™ = / da {t) (pW (f))> = £j dt (# (t))>, (47) 

= [da (t) (t) pw (t) = r°° dt^) (t) ^) x (t) ^-H 

J J — OO 

(48) 

for which the alternative explicit integral representations are shown to exist 
in the Appendix A. Without going into details of those calculations, we 

4 Hereafter we use the notation vn (A) for averaged density of states, Eq. (26). 
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only stress an extremely important role played by both the large-iV limit, 
Eq. (124), of the recurrence equation for associated orthogonal polynomials 
and the asymptotic expansion, Eq. (125), deduced from Eq. (124). We also 
notice that owing to the existence of these explicit representations, derived 
by using a large- N version of the recurrence equation (28), all further 
calculations became possible. 

The functions A$ and B^f 1 entering Eq. (46) can be expressed in terms 
of A { a N) and ri N) . Namely, bearing in mind the definitions given by Eqs. 
(34) and (35) we obtain for (A) = Ag? (A) + aA^ g (A), 



Ag? (A) = 2c Ar £4£A$_ (7) A 2 *- 2 , (49) 

k=l cr=l 

42(A) = 2 W /^(F<T»W) 2 . (») 
In the same way, the function B$> (A) = Bg? (A) + aSg (A) is given by 

<)(A) = 2c N ±d k j:r^i a) _ 1 X^-\ (51) 
k=l <t=1 



CiW = ^/^? (*)*(*)■ (52) 

When deriving these formulas we have used the fact of evenness of the 
measure da (t) jdt and of (pffl (i)) , as well as the expansion 

fk _ \k k 

— = V t k ~ m X m ~ 1 . (53) 

t — A -, 

m=l 

The "singular" components, A^2 g and B^^, can easily be determined. 
First, due to oddness of the integrand in Eq. (50), we have (A) = 0. 
Second, in order to find B^, we notice that the quantity 



In 



sing 



where n is not necessarily large, obeys the identity 7^ + 7^l\ = 1, which 
is a direct consequence of the recurrence equation (28). As far as 7^ = 
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72n-2 = • • • = 72°^ = 0) we conclude that 7^ = [1 — (—1)™] /2 and there- 



fore, 5^2 (A) = 2jffi /\. Hence we are led to the following representations, 



AP(X) = 2c^4EA^)^- 2 , (55) 

fe=l CT=1 

i#>(A) = 2c jV E4E r ^L)-i^- 1 +« \ j • (56) 

fe=l (7 = 1 

With these preliminarily calculations in hand we are ready to imple- 
ment the idea of recovering the Dyson density of states from the recurrence 
equation. In accordance with Eqs. (44) and (46) there are five contribu- 
tions to the integral in the r.h.s. of Eq. (44) corresponding to five terms in 
Eq. (46). Substituting Eqs. (55) and (56) into Eq. (46) and performing a 
formal integration with the help of Eqs. (47) and (48), these contributions 
are found to be 

Pl = 2c^£4EA^-^;li), (57) 

k=l a=l 

P2 = 2c%±d k ±A%Zy$ +s _ 1) , (58) 

k=l cr=l 

k=l <r=l 

k=l cr=l 

-ac N [l - (-lf]ri% (60) 

p fe-i 

P5 = 2CJVCJV-1 £ d k Yl r 2(fc-2)-i r 2(2f S )-l 
k=l a=\ 

+acw [i-(-ir]r<2, (6i) 

For large- N matrix models with a single spectrum support there are asymp- 
totic identities cat ~ cjv-i, A^ ~ A^ ^ and F^ ~ ^ (see Ap- 
pendix A for details) which simplify matters greatly. Collecting Eqs. (57) 
- (61) we come down to 

/-\-oo 
d\X 2s u N (A) -2ac N (-if 4^ 

-00 
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k=l cr=l 

Further double summation over indices k and a is performed by making 
use of the integral representations for and given by Eqs. (129) 
and (131) in Appendix A. Straightforward calculations lead to 

/+oo rT^N 
d\\ 2s v N (A) N ^°° / d\\ 2s Ps (A) , (63) 

with T>n = 2cat and 

-2L (64) 

For A = Patz with \z\ < 1 the term proportional to Xdv/dX vanishes iden- 
tically due to the principal value V of the integral over variable t, while the 
term proportional to a is a subleading one in the large-A?" limit. Then, in 
accordance with our concept Eq. (44), we end up with 

smooth (x) _ 2 / 2 2 y/2 f>» tdt dv/dt 

v n (A) [V N A ) VJ o (65) 

Equation (65) is exactly the Dyson density vd (A) of states with T>n being 
the end point of the eigenvalue support. We reconstructed the macroscopic 
level density Eq. (65) directly from the recurrence equation (28), alterna- 
tively to the traditional mean-field-theory derivation [1]. Notice that the 
spectrum end point T>n is the positive root to the integral equation 

ttN f v N tdt dv , . 

(66) 



2 Jo (v 2 N -t 2 ) 1/2 dt 
following from the normalization of the level density. 

4.2. GLOBAL CONNECTED "DENSITY-DENSITY" CORRELATION 
FUNCTION 

The same technology is applicable to the study of the smoothed connected 
"density-density" correlator pi N ^ . It is defined in terms of the scalar kernel 
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by Eq. (27), so that 

pf)(A,A') = _*2 (AjA , ) = __J^ 



(A)) (AO) + (A)) (# (AO) 

-2# (A) yjv-i (A) </>^ (AO^i(A')] • (67) 

We remind that here ipff (A) = exp{— V Q (A)} pj^ (A) are fictitious wave 
functions, A ^ A', and (. . denotes averaging over rapid oscillations mani- 
fested on the characteristic scale of the mean level spacing. The averaging 
in Eq. (67) can be done along the lines of the previous Subsection with 
two modifications. First, as far as A / A' we can run averaging over A 
and A' independently (this is already reflected in Eq. (67)). Second, we 

have to take into account the evenness of {^p^ (A)) and the oddness of 
There are two integrals 

r+oo 



A N) = |^°° ( iAA 2s e X p{-2V;(A)}(p^i(A)) 2 , (68) 
A N) = f +C °dXX 2s+1 exp{-2V a (X)}P i N a) (X)P i N a l 1 (X) (69) 

J — oo 

to be evaluated in the large- A" limit. With the help of Eqs. (47) and (48) one 

immediately recognizes them as the objects and rj,^^, respectively, 
calculated in Appendix A. By comparing of Eqs. (68) and (69) with Eqs. 
(129) and (131) we deduce that 



(„!?(A)) 2 = 1 -(v%-^y ,/2 , (70) 

vWmAm = (v% - aT" 2 , in) 

for |A| < V N . Substituting Eqs. (70) and (71) into Eq. (67) is a final step 
leading us to the smoothed "density-density" correlator 

AN) , x \/\ _ 1 V N ~ AA ; . . 

Pc l ' j " 2,Hx-xf {v %-x^ 2 (v%-x'^ 2 [ ] 

announced by Eq. (2) with (5 = 2. We stress, that it has been obtained here 
from the recurrence equation for orthogonal polynomials associated with 
the random matrix ensemble in question. 
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4.3. LOCAL EIGENVALUE CORRELATIONS BY SOLVING EFFECTIVE 
SCHRODINGER EQUATION 

4.3.1. Effective Schrddinger Equation 

Local eigenvalue correlations in the matrix ensemble Eq. (4) can be stud- 
ied by using the asymptotic version of the effective Schrodinger equation 
(36) obtained in Section 3. For the confinement potential V a introduced by 
Eq. (41) we obtain in the large-iV limit the following expressions for the 

functions F { n a) (A) and Q { n ] (A) (see Eqs. (37) and (38)), 



^ (A) 



Q { n ] (A) 



d_ 

dX 



(A), 



(73) 



i (a) 



(- 



d 



(a) 



a (—1) — a 
A 1 



(74) 



Here we have used the sum rule Eq. (39) to eliminate For confinement 
potential with a smooth regular part v, the second term in Eq. (74) is a 
subleading and hence it must be discarded 5 . Note, that it is rather inter- 
esting that the confinement potential V a does not appear in both equations 
above in an explicit way. It is even more exciting that in the considered ap- 
proximation the function Affi can be solely expressed through the Dyson 
density, Eq. (65). Indeed, taking into account the representation Eq. (132) 



and the fact that A$ (A) = AY^ (A), we are led to the asymptotic relation 



4? (A) = 



■kv d (A) 



(\/v N y 



1/2 ' 



(75) 



where, in accordance with Eq. (65), 



2 r v N 



tdt dv (V 2 N - A 2 \ 



t 2 - A 2 dt 



1/2 



v%-t 2 ) 



(76) 



This allows us to arrive at the following remarkable effective one-particle 
Schrodinger equation for the wave functions 



^ (A) = |A| Q P^ (A)exp{-t,(A)} 



(a) 



(77) 



5 This is not the case for the multicritical correlations near the origin A = 0. A detailed 
discussion of this important situation can be found in the very recent paper of Ref. [66] . 
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of fictitious non-interacting fermions in the large- N limit [48, 53] 



dX 2 



dX 



log 



■kv d (A) 



i - (x/v N y 



1/2 



A («) 



+ I^d(A) + ^ )^a/ ; (A) = 0. 



(78) 



Also, due to Eq. (29), one can verify that the wave functions y'iv-i (A) an d 
ipffi (A) °f two successive quantum states are connected by the relationship 



a O) 
dX 



TtV D (A) 



(A/Pjv) 5 



1/2 



(a) /n A ( a ) \ a (a) / x x 

¥>V-i (A) - pT^V (A)J + (-1) ^vV(A)- 



(79) 

Equations (78) and (79) serve as a general basis for the study of eigenvalue 
correlations in non-Gaussian random matrix ensembles in an arbitrary spec- 
tral range. 

It is instructive to analyze them in the particular case of GUE, where 
the Dyson density of states is the celebrated semicircle, 



; GUE 
D 



(A) = 7T- 1 (V% - X 



,y/2 



(80) 



with T>n = (2N) 1 / 2 . The square-root law for z/^ UE (A) immediately removes 

the first derivative d(pffi /dX in Eq. (78), providing us with the possibility 
to interpret the fictitious fermions as those confined by a quadratic poten- 
tial (a = 0). As far as the semicircle is a distinctive feature of density of 
states in GUE only, one will always obtain a first derivative in the effective 
Schrodinger equation for non-Gaussian unitary ensembles of random ma- 
trices. Therefore, fictitious non-interacting fermions associated with non- 
Gaussian ensembles of random matrices occur in a non-Hermitean quantum 
mechanics. 

An interesting property of these equations is that they do not contain 
the regular part of confinement potential explicitly, but only involve the 
Dyson density vd (analytically continued on the entire real axis) and the 
spectrum end point T>n- in contrast, the logarithmic singularity (that does 
not affect the Dyson density) introduces additional singular terms into Eqs. 

(78) and (79), changing significantly the behavior of the wave function tpffi 
near the origin A = 0. The influence of the singularity decreases rather 
rapidly outward from the origin. 

Structure of the effective Schrodinger equation leads us to the following 
statements [48] valid in the thermodynamic limit: 
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• Eigenvalue correlations are stable with respect to non-singular defor- 
mations of the confinement potential. 

• In the random matrix ensembles with well behaved confinement po- 
tential the knowledge of Dyson density (that is rather crude one-point 
characteristics coinciding with the real density of states only in the spec- 
trum bulk) is sufficient to determine the genuine density of states, as 
well as the n-point correlation function, everywhere. 

The latter conclusion is rather unexpected since it considerably reduces 
the knowledge required for computing n-point correlators. 

Effective Schrodinger equation obtained above enables us to examine 
in a unified way the local eigenvalue correlations in non-Gaussian ensem- 
bles with U (N) symmetry in different scaling limits. As we show below, 
it inevitably leads to the universal Bessel correlations in the origin scaling 
limit [47, 17, 53], to the universal sine correlations in the bulk scaling limit 
[36, 25, 38], and to the universal Airy correlations in the soft-edge scaling 
limit [53]. Corresponding scalar kernels are given by Eqs. (85), (87) and 
(94), respectively. 

4.3.2. Origin scaling limit and the universal Bessel law 
Origin scaling limit deals with the region of the spectrum close to A = 
where confinement potential displays the logarithmic singularity. In the 
vicinity of the origin the Dyson density can be taken as being approximately 
a constant, (0) = 1/A^r (0), where A^r (0) is the mean level spacing at 
the origin in the absence of the logarithmic deformation of potential v. 
Within the framework of this approximation, Eq. (78) reads 



Solution to this equation that remains finite at A = can be expressed by 
means of Bessel functions 



where a and b are constants to be determined later, and A (0) = A27V (0) ~ 
A2Ar+i (0). Inserting these solutions into Eq. (23) we find that the scalar 
kernel can be written down as 




(81) 




(82) 



(83) 
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-■4+1/2 ( ^y) 4-1/2 



(84) 



where the unknown factor c can be found from the requirement k!^ ^ (A, A) 
= 1/A(0). This immediately yields us the value c = — 7r/A(0). Defining 
now the scaled variable s = A s /A(0), we obtain that in the origin scal- 
ing limit the scalar kernel K OI i g (s, s') = hniN^oo K^2N ^s') takes the 
universal Bessel law, 

„ / a _ VT . 1/2 4+1/2 (7[£) 4-1/2 (tTsQ - 4-1/2 (^fi) 4+1/2 (^fiQ 
A orig ( s ) s )- 2 l SS J _ / 

(85) 

Equation (85) is valid for arbitrary a > —1/2, thus extending a recent proof 
[17] of universality of the Bessel kernel. 

4.3.3. Bulk scaling limit and the universal sine law 

Bulk scaling limit is associated with a spectrum range where the confine- 
ment potential is well behaved (that is far from the logarithmic singularity 
A = 0), and where the density of states can be taken as being approximately 
a constant on the scale of a few levels. In accordance with this definition 
one has 

Kbulk (s,s')= lim K olig (s,s'), (86) 

where s and s' should, nevertheless, remain far enough from the end point 
T>n of the spectrum support. 

Taking this limit in Eq. (86), we arrive at the universal sine law 

IT (S — S') 

deeply connected to the Wigner-Dyson level statistics [32]. 

4.3.4. Soft-edge scaling limit and the universal Airy law 

Soft-edge scaling limit is relevant to the tail of eigenvalue support where 
crossover occurs from a nonzero density of states to a vanishing one. It 
is known [50, 67] that by tuning coefficients dk which enter the regular 
part v of confinement potential (see Eq. (42)), one can obtain a macro- 
scopic (Dyson) density of states which possesses a singularity of the type 

vd (A) oc (1 — A 2 /V 2 N ) m+l l 2 with the multicritical index m = 0,2,4, etc. 
(Odd indices m are inconsistent with our choice that the leading coefficient 
dp, entering the regular component v (A) of confinement potential, be posi- 
tive in order to keep a convergence of integral for partition function in 
Eq. (4)). It was shown in Ref. [53] within the Shohat method that as long 
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as the multicriticality of the order m is reached, the eigenvalue correlations 
in the vicinity of the soft edge become universal, and are independent of 
the particular potential chosen. The order m of the multicriticality is the 
only parameter which governs spectral correlations in the soft-edge scaling 
limit. Here, however, we restrict ourselves to a general confinement poten- 
tial without tuning to the multicritical point, that corresponds to m = 0. 

Let us move the spectrum origin to its end point T>n, making the re- 
placement 



1 + 



(a/2) (2(ttP^ jv (1))- 1 



2/3" 



(88) 



,2/3 



that defines the soft-edge scaling limit provided a <C {V^IZn iX)Y' oc 
A 2 / 3 . It is straightforward to show from Eqs. (78) and (79) that in the 



vicinity of the end point T>n the function (p^ (s) 
the universal differential equation 



<Pn ( A s 



$N i s ) ~ si PN («) 

and that the following relation takes place, 



0, 



( I 



1/3 



<p' N (a) 



T>n) obeys 
(89) 

(90) 



Solution to Eq. (89) which decreases at s — > +oo (that is at far tails of 
the density of states) can be represented through the Airy function 



Ai( S ) = 3 



as follows 



s l/2 
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p N (s) = aAi (s) , 

with a being an unknown constant. Making use of Eq. (90), we obtain that 
in the vicinity of the soft edge the scalar kernel is 



Kn (A s , X s 



Ai (s) Af (a') - Ai (a') Ai' (a) 



s — a 



(93) 



where b is an unknown constant again. It can be found by fitting [51] 
the density of states Kn (A s , A s ), Eq. (93), to the Dyson density of states 
V D (A s ), Eq. (76), near the soft edge provided 1 <C s <C A 2 / 3 . This yields us 
the value b = c^ 1 (7tcjv71.jv(1)) 2//3 ■ Thus, we obtain that in the soft-edge scal- 
ing limit, Eq. (88), the scalar kernel Asoft {s, s') = limjv^oo \' S K^ (A s , A s ') 
satisfies the universal Airy law 



A so ft (s, a') 



Ai (a) Ai' (a') - Ai (a') Ai' (a) 



(94) 
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which does not depend on the details of the confinement potential. In fact, 
the Airy law is a particular case (m = 0) of more general multicritical corre- 
lations characterized by the index m of the multicriticality. For more details 
we refer the reader to Refs. [53, 48]. 

It follows from Eq. (94) that the density of states in the same scaling 
limit 

fsoft(s)= (^Ai(s)) 2 - S [Ai( S )] 2 (95) 

is also universal. The large-|s| behavior of zv so ft can be deduced from the 
known asymptotic expansions [68] of the Bessel functions, 

f N 1/2 1 rns f4|,3/2\ 

1 8^ eX P(-| s3/ > S^+OO. 

Note that the leading order behavior as s — > — oo is consistent with the 

1 /2 

\s\ 1 singularity of the bulk density of states. 
4.4. DISCUSSION 

Looking back at the formalism developed we should reiterate that the cru- 
cial point in the derivations above is the large- N limit, Eq. (124), of the 
recurrence equation for associated orthogonal polynomials. It was precisely 
this limit that led us to the important relation Eq. (75) and to the effec- 
tive Schrodinger equation in the form of Eq. (78) which is a nonuniversal 
in general. However, it takes locally universal forms in the spectrum bulk 
(where vn is approximately a constant on the scale of a few eigenlevels) , 
near the spectrum origin (where all the nontrivial information is contained 
in A~ 2 term in front of tpffi in Eq. (78)), and near the soft edge of the spec- 
trum (where universality shows up in the universal square-root singularity 
of I'd)- These three locally universal features of Eq. (78) have led us to the 
universal sine, Bessel and Airy kernels in corresponding scaling limits. 

We stress that Eq. (124) is, in fact, the leading-order-limit as N — > oo. 
How accurate is this approximation, and do situations exist where the next- 
order terms in the recurrence equation should be taken into account? A 
partial answer to this question was given in Ref. [66] whose authors, re- 
maining within the framework of the Shohat method, convincingly demon- 
strated that corrections to Eq. (124) are of importance in a problem of 
multicritical spectral correlations near the spectrum origin. The effective 
Schrodinger equation obtained there for two particular matrix ensembles 
with fine-tuned confinement potentials was shown not only to depend on 
the macroscopic spectral density vd (A) but, in addition, to contain con- 
tributions from subdominant terms in 1/N expansion for the recurrence 
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coefficients. It is important however, that in the situation in question the 
resulting differential equation contained the universal functions Fn and 
Qn involving certain universal combinations of recurrence coefficients and 
coupling constants responsible for the fine tuning of the multicritical con- 
finement potential. Having a complicated form, the effective Schrodinger 
equation could not be solved analytically, but, remarkably, the authors of 
Ref. [66] succeded in identifying a certain "mesoscopic" limit, in which the 
numerical solution of the exact differential equation and the analytical so- 
lution of the approximate differential equation obtained by making use of 
the relation Eq. (75) were shown to have quite similar qualitative features. 
With increasing of the order of the multicriticality near the spectrum origin, 
the approximate (analytical) and exact (numerical) solutions were shown 
to approach each other even quantitatively, demonstrating thus the poten- 
tialities of the Shohat method even in its simplest formulation presented 
above. 

5. Two Band Random Matrices 

5.1. MULTI-BAND SPECTRAL REGIMES 

Ensembles of large random matrices H generated by the joint distribu- 
tion function P[H], Eq. (1), may display phase transitions under non- 
monotonic deformation of the confinement potential V [H] . Different phases 
are characterized by topologically different arrangements of eigenvalues in 
random matrix spectra that may have multiple-band structure. Random 
matrices, whose spectra undergo phase transitions, appear in quantizing 
two-dimensional gravity [69, 70, 71], in the context of quantum chromody- 
namics [72, 73], as well as in some models of particles interacting in high 
dimensions [74] . Transition regimes realized in invariant random matrix en- 
sembles have implications for a certain class of Calogero-Sutherland-Moser 
models [75] . These matrix models may also be applicable to chaotic systems 
having a forbidden gap in the energy spectrum. 

It is convenient to parametrize the confinement potential V (A) entering 
Eq. (1) by a set of coupling constants {d} = {d\, ...,d p }, 



so that we may consider the phase transitions as occurring in {ci}-space. 
Because the confinement potential is an even function, the associated ran- 
dom matrix model possesses so-called Z2-symmetry. 

Variations of the coupling constants affect the Dyson density vd, that 
can be found by minimizing the free energy Fn = — log Zjy, Eq. (1), subject 




(97) 
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to a normalization constraint / vd (A) d\ = N, 

where T 7 indicates a principal value of the integral. When all dk are pos- 
itive, so that confinement potential is monotonic, the spectral density 
has a single-band support, A/b = 1. Non-monotonic deformation of the con- 
finement potential can be carried out by changing the signs of some of 
(k 7^ p). Such a continuous variation of coupling constants may lead, under 
certain conditions, to a discontinuous change of the topological structure 
of spectral density vp, when the eigenvalues {A} are arranged in A/b > 1 
"allowed" bands separated by "forbidden" gaps. 

The phase structure of Hermitean (/? = 2) one-matrix model Eq. (1) 
has been studied in a number of works [76, 77, 78, 79], where the simplest 
examples of non-monotonic quartic and sextic confinement potentials have 
been examined. It has been found that there are domains in the phase 
space of coupling constants where only a particular solution for i/£> exists, 
and it has a fixed number A/b of allowed bands. However, in some regions 
of the phase space, one can have more than one kind of solution of the 
saddle-point equation Eq. (98). In this situation, solutions with different 

number of bands Af^\ M^ 2 \ . . . are present simultaneously. When such an 

overlap appears, one of the solutions, say A/^ , has the lowest free energy 

(k) 

Fjy , and this solution is dominant, while the others are subdominant. 
Moreover, numerical calculations [78] showed that some special regimes 
exist in which the bulk spectral density obtained as a solution to the saddle- 
point equation Eq. (98) differs significantly from the genuine level density 
computed numerically within the framework of the orthogonal polynomial 
technique. It was then argued that such a genuine density of levels cannot 
be interpreted as a multi-band solution with an integer number of bands. 
A full understanding of this phenomenon is still absent. 

Recently, interest was renewed in multi-band regimes in invariant ran- 
dom matrix ensembles. An analysis based on a loop equation technique 
[20, 55] showed that fingerprints of phase transitions appear not only in the 
Dyson density but also in the (universal) wide-range eigenvalue correlators, 
which in the multi-band phases differ from those known in the single-band 
phase [19, 36, 23]. A renormalization group approach developed in Ref. [80] 
supported the results found in Refs. [20, 55] for the particular case of two 
allowed bands, referring a new type of universal wide-range eigenlevel cor- 
relators to an additional attractive fixed point of a renormalization group 
transformation. 

As it was already stressed in the Introduction, the method of loop equa- 
tions [20, 55], used for a treatment of non-Gaussian, unitary invariant, 
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random matrix ensembles fallen in a multi-band phase, is only suitable for 
computing the global characteristics of spectrum. Therefore, an appropriate 
approach is needed capable of analyzing local characteristics of spectrum 
(manifested on the scale of a few eigenlevels) . A possibility to probe the 
local properties of eigenspectrum is offered by the method of orthogonal 
polynomials. A step in this direction was taken in the paper [81], where 
an ansatz was proposed for large-V" asymptotes of orthogonal polynomi- 
als associated with a random matrix ensemble having two allowed bands 
in its spectrum. Because the asymptotic formula proposed there is of the 
Plancherel-Rotach type [60], it is only applicable for studying eigenvalue 
correlations in the spectrum bulk and cannot be used for studying local 
correlations in an arbitrary spectrum range (for example, near the edges of 
two-band eigenvalue support). 

Below we demonstrate that the Shohat method needs minimal modifi- 
cations to allow a unified treatment of eigenlevel correlations in the unitary 
invariant U (N) matrix model ((3 = 2) with a forbidden gap. In particular, 
we will be able to study both the fine structure of local characteristics of the 
spectrum in different scaling limits and smoothed global spectral correla- 
tions. As is the single-band phase, the treatment presented below is based 
on the direct reconstruction of spectral correlations from the recurrence 
equation for the corresponding orthogonal polynomials. 

5.2. EFFECTIVE SCHRODINGER EQUATION IN THE TWO-BAND 
PHASE AND LOCAL EIGENVALUE CORRELATIONS 

Let us consider the situation when the confinement potential has two deep 
wells leading to the Dyson density supported on two disjoint intervals lo- 
cated symmetrically about the origin, V N < |A| < V N . In this situation, 
the recurrence coefficients c n entering Eq. (28) are known to be double- 
valued functions of the number n [69, 79]. This means that for n = N » 1 
and in contrast with a single-band phase, one must distinguish between 
coefficients Cjv±2g ~ Qv and coefficients CN-\±2q ~ cjv-i, belonging to two 
different, smooth (in index) sub-sequences; here, integer q ~ O (N°). Bear- 
ing this in mind, the large-V version of recurrence equation Eq. (28) can 
be rewritten as 

A 2 - (4 + c 2 v _ 1 )] P N (A) = Civciv-i [Pn-i (A) + P N +i (A)] , (99) 

whence the two analogues of the asymptotic expansion Eq. (125) can be 
obtained. They are given by Eqs. (135) and (136) of Appendix B. We no- 
tice that this is the only point crucial for extending the Shohat method 
to the double-well matrix models considered in this Section. These new 
expansions make it possible to compute the required functions and On 
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entering the differential equation Eq. (36) for fictitious wave functions in 
the limit N » 1. 

In accordance with the general framework of the Shohat method, we 
have to compute two functions (compare with Eqs. (55) and (56)) 

v k 



A N {\) = 2c N J2 d kJ2 A 
B N (A) 



fc=i 
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(100) 
(101) 



involving the objects and r^j^, for which there exist the useful in- 
tegral representations given by Eqs. (140) and (148) in Appendix B. Sub- 
stituting them into Eqs. (100) and (101) one is able to perform the double 
summation over indices k and a. Omiting details of straightforward calcu- 
lations we present the final answer given by the formulas 
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(103) 

Note, that along with a different (compared to the single-band phase) func- 
tional form of the functions Am (A) and Bn (A), these functions are, in fact, 
double-valued in index N, and behave in a different way for odd and even 
N. This is a direct consequence of the "period-two" behavior [69, 79] of 
the recurrence coefficients c n . 

Having obtained the explicit expressions for functions An and Bjy, it 
is easy to verify that coefficients T n (A) and Q n (A) entering the differential 
equation Eq. (36) for the fictitious wave function ip n (A) may be expressed 

in terms of the Dyson density vffl in the two-cut phase supported on two 
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when N S> 1. This formula can be obtained either via the procedure of 
the Sec. 4.1 or within the mean-field approach, Eq. (98). Here and 
are the end points of the eigenvalue support that obey the two integral 
equations 
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derived in Appendix C. 

By making use of the Eqs. (37), (38) and (102) - (104) we obtain in the 
large- N limit 
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so that for N » 1 the effective Schrodinger equation in the two-cut phase 
reads [56] 
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(109) 



As tends to zero, we reproduce the equation (78) with a = valid in 
the single-band regime. 

Local eigenvalue correlations in the spectra of two-band random matri- 
ces are completely determined by the Dyson density of states entering the 
effective Schrodinger equation Eq. (109). 

(i) In the spectrum bulk, the Dyson density is a well behaved function 
that can be taken approximately as being a constant on the scale of a few 
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eigenlevels. Then, in the vicinity of some Ao that is chosen to be far enough 
from the spectrum end points ±X>^, Eq. (109) takes the form 

^g^ + [vr/A(A )] 2 ^(A) = 0, (110) 

with A (Ao) = 1/^d^ (Ao) being the mean level spacing in the vicinity of 
Ao- Clearly, the universal sine law, Eq. (87), for the two-point kernel follows 
immediately. 

(ii) To study the eigenvalue correlations near the end points of an eigen- 
value support we notice that in the absence of the fine tunning of confine- 
ment potential, the Dyson density has a universal square-root singularity 

in the vicinity of |A| = V%, that is (A) oc ^1 - (x/V^^j . We 

then readily recover the universal Airy correlations, Eq. (94), previously 
found in the soft -edge scaling limit for U (N) invariant matrix model in 
the single-band phase. 

5.3. GLOBAL CONNECTED "DENSITY-DENSITY" CORRELATOR 

Let us turn to the study of the smoothed connected "density-density" cor- 
relator that is expressed in terms of the scalar kernel as follows (see Eq. 
(27)), 



P2P(A,V) = -^^{^(A)^ 1 (A') + ^(A')^-i(A) 

-2<p N (A) ip N -! (A) ip N (A') (A')} • (HI) 

Here A 7^ A'. Equation (111) contains (before averaging) rapid oscillations 
on the scale of the mean level spacing. These oscillations are due to presence 
of oscillating wave functions (fN and ifiN-i- 

To average over the rapid oscillations, we integrate, over the entire real 
axis, rapidly varying wave functions in Eq. (Ill) multiplied by an arbitrary, 
smooth, slowly varying function, which without any loss of generality can 
be choosen to be X 2s for ip 2 N (A) and A 2s+1 for yjv (A) yjv-i (A) (s is an 
arbitrary positive integer, s > 0). Consider, first, the integral 

l[ N) = f + °°dXX^ N (X)=A^. (112) 
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With the help of Eq. (140), and bearing in mind that (p% (A) is an even 
function, we conclude that 

IAI X 2s d\ 
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whence, in the large- N limit, 
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with 6 being a step function. The same procedure should be carried out 
with the expression ip^ (A)</?jv-i (A) in Eq. (111). Since this construction 
is an odd function of A, we have to consider the integral 
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With the help of Eq. (148), and exploiting the oddness of ip n (A) fN-i (A), 
we rewrite Eq. (116) in the form 
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Combining Eqs. (Ill), (114), (118) and (147), we finally arrive at the fol- 
lowing formula for smoothed "density-density" correlator [56] 
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where 
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It is seen from Eq. (120) that smoothed "density-density" correlator in 
the two-band phase differs from that in the single-band phase, Eq. (72). 
However, it is still universal in the sense that the information of the distribu- 
tion Eq. (1) is encoded into the "density-density" correlator only through 
the end points of the eigenvalue support. The striking parity effect 
in the new universal function Eq. (120), that is the sharp dependence of 
correlations on the oddness/evenness of the dimension N of the random 
matrices, is the main qualitative difference as compared to the global core- 
lations in random matrices fallen in the single-band phase. This effect is 
most pronounced in the case of unbounded spectrum. The origin of this 
unusual large- N behavior will be discussed later on. 

Finally, let us speculate about the universal correlator Eq. (120) in 
the limit of unbounded spectrum, — ► oo, with a gap. Inasmuch as it 
describes correlations between the eigenlevels which are repelled from each 
other in accordance with the logarithmic law, that is known to be realized 
[82, 83] in weakly disordered systems on the energy scale |A — A'| <C E c (E c 
is the Thouless energy), we may conjecture that the corresponding limiting 
universal expression 

lim pg?(A,A') = - y (AA 2 9(|A|-A)9(|A1-A) 

2ir 2 (A - A') 

AA' - A 2 
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reflects the universal properties of real chaotic systems with a forbidden 
gap A = and broken time reversal symmetry, provided |A — A'| <C E c . 
In two limiting situations (i) of gapless spectrum, A = 0, and (ii) far from 
the gap, |A| , |A'| A, the correlator Eq. (121) coincides with that known 
in the Random Matrix Theory of gapless ensembles [36, 23] and derived in 
Ref. [82] within the framework of diagrammatic technique for spectrum of 
electron in a random impurity potential. 



5.4. DISCUSSION 

In this Section we have demonstrated how the Shohat method should be 
transformed in order to study both global and local spectral characteris- 
tics of U (iV) invariant ensembles of large random matrices possessing Z2- 
symmetry, and deformed in such a way that their spectra contain a forbid- 
den gap. We proved that in the pure two-band phase, the local eigenvalue 
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correlations are insensitive to this deformation both in the bulk and soft- 
edge scaling limits. In contrast, global smoothed eigenvalue correlations in 
the two-band phase differ drastically from those in the single-band phase, 
and generically satisfy a universal law, Eq. (120), which is unusually sen- 
sitive to the oddness/evenness of the random matrix dimension provided 
the spectrum support is bounded. On the formal level, this sensitivity is 
a direct consequence of the "period-two" behavior [69, 79] of the recur- 
rence coefficients c n that is characteristic of two-band phase of reduced 
Hermitean matrix model. To see this, consider the simplest connected cor- 
relator (TrHTrH) c that can be exactly represented in terms of recurrence 
coefficients for any n, 

(Tr HTr H) c = c\. (122) 

Since in the two-band phase double-valued function of index n, 

alternating between two different functions as n goes from odd to even, the 
large- N limit of the correlator (TrHTrH) c strongly depends on whether 
infinity is approached through odd or even N. Then, an implementation of 
a double- valued behavior of c n into the higher order correlators of the form 
^TrH fc TrH'^> contributing to the connected "density-density" correlator 

gives rise to the universal expression Eq. (120), which is valid for the two- 
band random matrix model with pure Z2-symmetry. 

Let us, however, point out that no such sensitivity has been detected in 
a number of previous studies [20, 55] exploiting a loop equation technique. 
One possible explanation comes from the following reasons. In the method 
of loop equations, used for a treatment of non-Gaussian random matrix 
ensembles fallen in a multi-band phase, one has to keep the most general 
(non-symmetric) confinement potential V (A) = Y^k=i dk^ k I 'k until very 
end of the calculations, and to take the thermodynamic limit N — > oo prior 
to any others. Therefore, Z2-symmetry in this calculational scheme can 
only be implemented by restoring Z2-symmetry at the final stage of the 
calculations, setting all the extra coupling constants c^fc+i to zero. Doing so, 
one arrives at the results reported in Refs. [20, 55]. From this point of view, 
our treatment corresponds to the opposite sequence of thermodynamic and 
Z2-symmetry limits, since we have considered the random matrix model 
that possesses Z2-symmetry from the beginning. Qualitatively different 
large-iV behavior of the smoothed connected "density-density" correlator, 
Eq. (120), and of the smoothed connected two-point Green's function given 
by Eq. (15) of Ref. [55] provides a direct evidence that the order of ther- 
modynamic and Z2-symmetry limits is indeed important when studying 
global spectral characteristics of multi-band random matrices. 

The parity effect manifested in global spectral correlators of double- 
well matrix models was the focus of the discussion in the recent study [84] . 
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The authors of Ref. [84] noted that, contrary to the method of orthogonal 
polynomials, the standard large-A limit techniques of analyzing matrix 
models like the loop equation method [19, 20] and the renormalization 
group approach [85] assume a smooth behavior with respect to N in the 
thermodynamic limit. The result Eq. (120) obtained by the authors of Ref. 
[84] in a different way [81] led them to conclusion that these methods need 
to be revisited when one deals with matrix models posessing eigenvalue 
gaps. 

6. Conclusion 

In this review we presented a formalism for statistical description of spectra 
of U (A) invariant ensembles of large random matrices. It lies within the 
general framework of the orthogonal polynomials' technique, and consists 
of the direct reconstruction of spectral densities and spectral correlations 
from the recurrence equation for orthogonal polynomials associated with 
a given random matrix ensemble. We have demonstrated the potentialities 
of this method, considering in a unified way both global and local spectral 
characteristics in matrix models with and without an eigenvalue gap. Al- 
though we directed our main attention to the most known bulk, origin and 
soft-edge scaling limits characterized by the universal sine, Bessel and Airy 
kernels, respectively, there are examples in the recent literature signaling 
about applicability of the described formalism to more refined situations - 
such as multicritical correlations near the soft edge of the spectrum support 
[53] and near the spectrum origin [66]. 

Attaching special significance to the study of the large- N limit of the 
recurrence equation for associated orthogonal polynomials, this method 
turns the recurrence equation into a kind of laboratory allowing the con- 
struction of matrix models with nonstandard properties - for example with 
eigenvalue gaps - by guessing a particular ansatz for the behavior of the 
recurrence coefficients cn in the thermodynamic limit. Just this feature of 
the formalism presented forces us, finally, to mention a crucial difference 
between the random matrix ensembles with strong level confinement con- 
sidered here and the random matrix ensembles with extremely soft level 
confinement. While the former ensembles (with confinement potentials of 
the Freud and Erdos type [38]) are characterized by a powerlike large- A 
limit of recurrence coefficients, cjv oc N p (p > 0), the latter (representing 
a class of g-deformed potentials [86, 87]) exhibit a qualitatively different, 
exponential rate of growth, cjy oc q N (q > 1). This results in a different 
large- N limit of the recurrence equation that will not already be as simple 
as stated in Eq. (124), and the emergence of different nontrivial classes of 
spectral statistics is inevitable. We consider a treatment of g-deformed ran- 
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dom matrix ensembles as a challenge to the Shohat method which, going 
back to 1930, had to wait so long to find its application in Random Matrix 
Theory. 
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A. Integral Representation of and r^L: Single Band Phase 

Consider the integral 

A™ = |da(f)(pW(t))> (123) 

with integer a > 0. In the large- A limit an alternative explicit integral 
representation can be found for A^P . This is achieved by making use of 
the large- A version of the recurrence equation (28). It is known that in the 
single-band phase of the matrix model the recurrence coefficients approach 
[69, 79] a smooth (in index A) single-valued function, so that c^+q ~ cn 
for q being of order O (A ). Within this approximation one obtains from 
Eq. (28) 

XP™ (A) = c N (P^ (A) + pW (A)) , (124) 
whence it follows that 

X m P^ (A) = ]T ( 7 ) P { N %_ m (A) , m > 1. (125) 



3=0 



J 



The identity (125) can be proven by the mathematical induction. The ad- 
vantage of the asymptotic expansion Eq. (125) is that being substituted 
(for m = 2a) into Eq. (123), it immediately allows us to explicitly perform 
the integration due to the orthogonality property Eq. (17). This yields 

a£° =<#!;( 2 ;H- d26) 



3=0 



with 6j being the Kronecker symbol. To evaluate the sum in Eq. (126) we 
make use of the integral representation of the Kronecker symbol, 
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We then find that 

A£° = # Re fj | E ( 2 ; ) exp {2i9 (j - a)} , (128) 

and hence, after summation over j and some rearrangements, 

rV» dtt 2 ° 

Here T>n = 2cjv. 
The integral 



2 [° N dtt 



a.™ 



4% =Jda (t) p£> (t) (t) t 2 ^ (130) 

with integer a > is computable in the same manner, with an answer given 
by 



W^Jo (V 2 N -t 2 ) l/2 ' 



(131) 



Notice that due to the asymptotic property c^+q ~ cjv mentioned above, 

a£^ « A£P and C « for g ~ O (JV°). 

Finally, we demonstrate the usefulness of the formulas (129) and (131) 

by finding the explicit expressions for the functions A^g (A) and Brej? (A), 
defined by Eqs. (49) and (51). Substitution of Eq. (129) into Eq. (49) fol- 
lowed by summation over a yields 

A (N) m ™ N f>« dt * t 2k - \ 2k 



77 (V 2 N -t 2 ) 1/2 t 2 - x2 ' 

Analogously, we obtain 

rW m 2 trff Acfo/rft - tdv/dX 

rcg 1 j vr 7 7o {V%-t 2 ) l/2 t 2 -^ 

2 A (it tdv/dt dv 

~ ~ V Jo (Tft N - t*) 1/2 * ~ >? ~ (133) 

One can convince himself that Eqs. (132) and (133) obey the sum rule Eq. 
(39) for the confinement potential V a , Eq. (41), and for the functions Aff 
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and B]£, given by Eqs. (55) and (56). Notice that this is the recurrence 
equation (28) which enabled us to obtain closed analytic expressions (132) 
and (133) relating the functions A^g and B^g to the regular part v of the 
confinement potential. 

B. Integral Representation of and r^]^: Two Band Phase 

Consider the integral 

a£°= J da{t)P 2 N {t)t 2 ° (134) 

with integer a > 0. It follows from Eq. (99) that in the two-band phase the 
following asymptotic identities exist, 




(135) 



and 

^iww - (*+*_)-£( J) 

X [cjV-iPjV+4j"-2fc+l W + CjV-FW+4j-2fc-l (A)] (136) 

with integer m > 0. Both Eqs. (135) and (136) can be proven by the 

mathematical induction. Making use of Eq. (135) we rewrite A 2 ^ in the 
form 




(137) 



Orthogonality of P n allows us to integrate over the measure da, thus sim- 
plifying Eq. (137) to 

^-(*^)'£(:)(^)s(5)*- (138) 
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Substituting the integral representation (127) for the Kronecker symbol, 
and performing the double summation over indices j and k, we obtain 



, (N) f 2w dO f o o „ \c 

A 2a = J o 2^ [ C N + C N-1 + 2c N C N -i COS 9 J 



(139) 



Introducing a new integration variable t 2 = c 2 N + c^ r _ 1 + 2ctvCat_ 1 cos 9, we 
derive the final formula 



A 



(AT) 
2a 



U N 



t 2 ° +1 dt 



1/2 



1/2 



with 

The integral 



T^n = \°n ± cn- 



4% = J da (t)P N (t)P N ^(t)t 



2a+l 



(140) 



(141) 



(142) 



with integer a > is evaluated in the same way. Making use of expansion 
Eq. (136), we rewrite Eq. (142) in the form that allows us to perform the 
integration over the measure da, 



2cr+l 



1 f ^ 

I, (c% + <?n-i) / da (t) Pa,.! (t) £ 
J k=0 



C 2 N + C Ar _ 1 



3=0 



[c N -lPN+4j-2k+l (*) + C N PN+4j-2k~l (*)] • (143) 



After integration, we get 

r L+i = 2 + c n-i) E ( 



fc=0 



CATCAr_i 



Cjy + C A r_ 1 



xE 

j=0 



CiV-l^+l + CN$2j 



(144) 



The double summation in Eq. (144) can be performed by using the integral 
representation for the Kronecker symbol given by Eq. (127), 



r 



(AO 

2(7+1 



1 f 2w d9 / \<r 

2 Jo 2tt \ n + C ^ _1 + 2cnCn ' 1 cos °) t CAr + CAr - 1 cos ^ ■ 



(145) 
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Introducing a new integration variable t 2 = c 2 N + cj^_i + 2cncn~i cos 6, we 
get 



r (JV) _ J_ f 



v+ 



t 2 ° +1 dt 



1/2 



1/2 



* 2 + C N ~ C N-1 



(146) 

Notice that because P_i (A) = 0, it follows from Eq. (28) that Co = 0, and 
as a consequence, an even branch C2N always lies lower than an odd branch 
C2JV±i> so that C2N < C2jv±i- Then, we may conclude from Eq. (141) that 



CAT 



V+-(-lfy- 



(147) 



and, as a consequence, 



1 2(7+1 — 



1 



t 2 °+ l dt 



1/2 



K) -* 2 



1/2 



(148) 



C. Soft Edges in the Two Band Phase 

To find the equations determining the end points T>^ where the Dyson 
spectral density goes to zero, we start with the following formula from the 
theory of orthogonal polynomials [63] 



da (t) — P n (t) P n _i (t) , 



(149) 



also known as a "string equation". Let us use expansion Eq. (136) to eval- 
uate the integral entering Eq. (149) in the limit n = N S> 1. It is easy to 
see that 

p r p 
N = 2c N Y / dkj da (t) P N (t) P N _ X (t) t 2 ^ 1 = 2c N £ d k T { £l v (150) 



fc=i 



where ^^-i 1S gi ven by Eq. (146). Then, we immediately obtain the rela- 
tionship 



N 



2 fK 



2 r 1 

TT Jv 



dt 



dV 



1/2 



{K) -t 2 *-{v- N ) 



V2 dt 



.2,2 2 



(151) 
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This result, rewritten for n = N — 1, yields in the large- N limit, 



N = — 

vr Jx>- 



dt 



dV 



1/2 r 



K) -t 2 ' t 2 -(v- N ) 



V2 dt 



t 2 + c^.]^ — c% 



(152) 

Equations (151) and (152) bring us two integral equations whose solutions 
determine the end points V^, 



f 

Jv~ 



v+ 



t 2 dt 



1/2 r 



V2 dt 



dV ttN 



(153) 



and 



vi 



dt 



dV 



v~ 



1/2 



{K) -t 2 t*-(v- N ) 



V2 dt 



0. 



(154) 



As Djy — > 0, Eq. (153) coincides with the integral equation (66) for a 
single-band phase. In the same limit, Eq. (154) becomes equivalent to the 
assertion U£> (0) = 0, with being the spectral density in a single-band 
phase. This corresponds to the point of merging of two eigenvalue cuts. 
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